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25  Abstract 


26  The  noctilucent  clouds  (NLC)  are  high-altitude  bright  cloud  formations  visible  under 

27  certain  conditions  from  high-latitude  places  during  the  summer  months.  Even  if  the  exact 

28  nature  of  these  clouds  still  remains  a  mystery,  they  are  an  efficient  tracer  of  the  dynamic 

29  processes  at  their  level,  particularly  the  gravity  waves  propagating  from  the  stratosphere 

30  through  the  mesopause  layer.  In  this  paper,  we  describe  a  technique  developed  to  analyze 

31  the  structures  visible  in  the  NLC  images  taken  every  summer  night  since  2004  from 

32  Stockholm,  Sweden  (59.4°N).  The  parameters  of  30  short-period  gravity  wave  events 

33  have  been  measured  and  compared  with  older  datasets  obtained  mostly  from  low  and 

34  mid-latitude  sites,  using  airglow  imaging  techniques.  The  horizontal  wavelengths  are  in 

35  good  agreement  with  previous  results  while  the  observed  horizontal  phase  speeds  exhibit 

36  smaller  values  than  for  other  sites.  The  directionality  of  the  waves  presents  strong 

37  poleward  preference,  traditionally  observed  during  the  summer  season.  This  anisotropy 

38  and  the  difference  in  the  phase  speed  distribution  cannot  be  explained  by  the  filtering  due 

39  to  the  background  wind  field  but  more  probably  by  the  position  of  the  gravity  waves 

40  sources,  located  to  the  south  of  the  observation  site. 

41 

42  1.  Introduction 

43  Noctilucent  clouds  (NLCs)  are  thought  to  be  composed  of  small  ice  particles  but  their 

44  process  of  formation  is  still  a  hot  research  topic.  They  occur  in  the  upper  mesosphere, 

45  between  80  and  85  km,  with  a  mean  altitude  of  82.9  km  (Gadsden  and  Schroder,  1989), 

46  and  are  most  frequently  observed  from  the  ground  between  50  and  65°  of  latitude.  Their 

47  appearance  is  very  similar  to  cirrostratus  clouds,  which  are  also  composed  of  ice  but 

48  develop  in  the  troposphere;  however  their  color  is  generally  silvery  blue  (Fig.  1).  Each 


49 


NLC  may  be  composed  of  different  forms,  classified  in  four  major:  veil  (Type  I),  bands 

50  (Type  II),  billows  or  waves  (Type  III)  and  whirls  (Type  IV).  These  four  major  types  are 

51  then  separated  in  several  subdivisions  (Gadsden  and  Parviainen,  1995).  NLCs  were  first 

52  reported  in  1884  (Leslie,  1884;  Jesse,  1889),  but  their  occurrence  and  brightness  are 

53  claimed  to  be  increasing,  especially  during  the  last  30  years  (Gadsden,  2002; 

54  Klostermeyer,  2002;  Deland  et  ah,  2003].  Furthermore,  sightings  from  unusual  mid- 

55  latitude  sites  have  more  often  been  reported  (Taylor  et  ah,  2002;  Wickwar  et  ah,  2002; 

56  Nielsen  et  ah,  this  issue).  Following  these  observations,  it  has  been  suggested  that  the 

57  occurrence  of  NLCs  may  be  a  very  sensitive  indicator  of  the  changes  in  the  atmosphere 

58  composition  and  dynamics,  and  possibly  of  the  Earth  global  change. 

59  Besides  the  study  of  their  own  properties  (Witt,  1962;  Fogle  and  Haurwitz,  1966), 

60  NLCs  are  also  an  interesting  tool  for  analyzing  the  dynamics  of  the  upper  atmosphere.  In 

61  fact,  they  work  as  a  tracer  of  the  coherent  structures  forming  or  propagating  at  their  level 

62  (e.g.,  Fritts  et  ah,  1993;  Chandran  et  ah,  2009).  Optical  studies  of  similar  phenomena 

63  have  already  been  done  using  the  mesopause  region  airglow  emissions  as  markers  of  the 

64  waves  (e.g.,  Taylor  et  ah,  1997;  Pautet  et  ah,  2005;  Nielsen  et  ah,  2009).  Since  the  NLC 

65  layer  is  significantly  thinner  than  a  typical  airglow  layer  (2-3  km  compared  to  8-9  km) 

66  (She  and  Lowe,  1998),  the  observed  wave  signatures  should  appear  much  more  clearly  in 

67  NLC  images  and  thus  could  be  studied  in  much  better  detail.  Also,  the  airglow  imaging 

68  technique  requires  conditions  that  are  not  always  encountered  at  high-latitudes  because  of 

69  the  persistent  twilight  in  the  summer  months  and  the  possible  auroral  contamination  of 

70  these  faint  emissions  throughout  the  year.  Hence,  imaging  NLCs  provide  a  method  to 

71  perform  summer-time  observations,  and  in  conjunction  with  airglow  observations, 
complete  a  full  seasonal  study  of  wave  structures. 
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73  Several  types  of  structures,  with  different  origins,  may  be  analyzed  using  the  patterns 

74  visible  in  the  NLC  images  (Haurwitz,  1961),  like  the  short-period  (<lh)  gravity  waves  (or 

75  bands)  and  the  ripples  (or  billows).  The  background  wind  might  also  be  determined 

76  (Fogle  and  Haurwitz,  1966,  Fritts  et  ah,  1993). 

77  The  small-scale  waves,  termed  ripples  or  billows,  are  short-lived  (typically  few  tens 

78  of  minutes)  structures  due  to  localized  regions  of  convective  or  shear  instabilities, 

79  generated  in-situ  at  mesospheric  altitudes  (Taylor  and  Hapgood,  1990;  Fritts  et  ah,  1993; 

80  Fritts  et  ah,  1994).  They  may  also  be  the  result  of  the  breaking  of  a  freely  propagating 

81  gravity  wave  (Horinouchi  et  ah,  2002). 

82  Gravity  waves  are  the  oscillations  of  air  parcels  by  the  lifting  force  of  buoyancy  and 

83  the  restoring  force  of  gravity  (Hines,  1960).  They  can  be  generated  by  several  sources:  in 

84  the  troposphere  by  thunderstorm  updrafts,  frontal  systems  or  airflow  over  mountains,  but 

85  also  by  interaction  with  the  polar  jet  stream,  geostrophic  adjustments,  wave-wave 

86  interactions  or  wave  breaking.  These  waves  propagate  vertically  and  horizontally  and 

87  may  be  observed  during  several  hours.  Their  horizontal  wavelength  may  vary  from  10  to 

88  more  than  100  kilometers.  They  play  an  essential  role  in  affecting  the  dynamics  and 

89  structures  of  the  middle  and  upper  atmosphere  and  participate  in  driving  the  thermal 

90  structure  and  setting  up  large-scale  meridional  circulation  in  the  mesosphere  and  lower 

91  thermosphere  (MLT)  region  (altitude  -60-100  km).  Freely  propagating  gravity  waves 

92  actively  transport  energy  and  momentum  from  the  troposphere  into  the  MLT  region, 

93  where  they  deposit  their  energy  and  transfer  their  momentum  to  the  mean  flow  when 

94  breaking  in  the  more  rarefied  air  (Lindzen,  1981;  Fritts  and  Alexander,  2003).  They  also 

95  have  direct  effects  on  the  noctilucent  clouds:  horizontal  variations  in  cloud  particle 
concentration,  decrease  of  ice  crystal  size  and  average  cloud  albedo  for  short-period  (less 
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97  than  a  few  hours)  gravity  waves,  inhibition  of  the  NLC  formation  due  to  temperature 

98  perturbations  (Jensen  and  Thomas,  1994;  Rapp  et  ah,  2002). 

99  When  seen  from  the  ground,  all  these  waves  are  distorted  by  the  geometry  of  the 

100  observation.  The  spherieal  shape  of  the  atmospherie  layer  where  NLCs  reside  and  the 

101  refraction  due  to  propagation  through  the  atmospheric  medium,  modify  the  shape,  the 

102  size  and  the  observed  speed  of  the  waves.  All  these  parameters  are  necessary  to 

103  understand  the  gravity  waves  properties,  thus,  it  is  important  to  work  with 

104  geographically-oriented,  linear-scaled  images  to  permit  the  accurate  measurement  of 

105  these  values. 

106  In  this  paper,  we  will  describe  a  technique  developed  to  analyze  short-period  gravity 

107  wave  structures  observed  in  NLC  images.  Following  a  detailed  description  of  the 

108  problem  geometry,  the  original  images  will  be  projected  onto  a  linear  geographic  grid  in 

109  a  process  known  as  “unwarping”.  These  new  images  will  be  similar  to  satellite -type 

110  pictures:  the  NLC  layer  as  seen  from  the  top,  looking  downward,  oriented  due  to  the 

111  geographic  north,  and  represented  with  a  linear  scale.  The  second  part  of  this  paper 

112  presents  an  analysis  of  the  short-period  gravity  waves  observed  during  5  NLC  seasons, 

113  between  2004  and  2008,  from  the  Stockholm  University  campus  (59.4°N).  Their 

114  characteristics  are  compared  with  previous  low  and  mid-latitude  and  recent  high-latitude 

115  measurements  obtained  using  airglow  data. 

116 

117  2,  Observations 

118  NLCs  may  be  seen  at  mid  and  high-latitudes  due  to  the  scattering  of  sunlight  from  a 

119  very  thin  layer  of  tiny  ice  particles  (~50  nm)  which  can  form  and  grow  in  the  very  cold 

120  conditions  (-'120-140  K)  of  the  high-latitude  mesopause  region  during  the  summer 


121  months.  In  this  specific  latitudinal  zone,  when  the  sun  lies  between  6  and  12°  below  the 

122  horizon  (nautical  twilight),  an  observer  staying  in  the  darkness  will  be  able  to  see  NLCs 

123  which  are  still  in  the  direct  sunlight  because  to  their  high  altitude.  Under  these  conditions, 

124  the  sunlight  scattering  from  the  atmosphere  occurs  only  at  higher  altitudes  (above  at  least 

125  40-50  km)  and  the  sky  is  dark  enough  to  observe  NLCs  that  are  seen  within  the  twilight 

126  arch  as  it  moves  in  azimuth  with  the  sun  (Fogle  and  Haurwitz,  1966). 

127  Since  summer  2004,  pictures  of  NLCs  have  been  taken  from  the  top  floor  window  of 

128  the  Arrhenius  Laboratory  at  the  university  campus  in  Frescati,  Stockholm,  Sweden 

129  (59.4°N,  18.1°E).  A  digital  camera  (Canon  PowerShot  G5),  operating  in  interval  time 

130  shooting  mode,  takes  every  night  hundreds  of  pictures  of  the  twilight  sky  at  the  rate  of  1 

131  (2004  and  2005)  or  2  (since  2006)  pictures  per  minute,  with  exposure  times  varying  from 

132  0.5  to  6s,  depending  on  the  darkness  of  the  sky.  The  images  are  stored  on  a  computer  in  a 

133  digital  format.  Fig.  1  shows  an  example  of  NLCs  photographed  during  the  night  of  16-17 

134  July,  2005  at  22:30  UT.  Many  structures  with  various  shapes  and  sizes  are  clearly  visible 

135  in  the  twilight  arch.  At  the  upper  edge  of  the  image,  above  the  twilight  arch,  the  NLCs  are 

136  not  visible  as  they  are  not  illuminated  by  the  sunlight,  and  in  the  lower  part  of  the  image, 

137  below  the  twilight  arch,  because  the  sky  is  too  bright  as  the  sunlight  reaches  lower 

138  atmospheric  layers.  It  is  possible  to  see  short-period  gravity  waves  in  the  center  right  half 

139  of  the  image,  at  mid  height,  but  also  ripples  structures  at  the  same  elevation,  close  to  the 

140  right  edge.  Although  other  structures  are  present  in  the  image  we  will  restrict  our  analysis 

141  to  these  two  groups  of  structures.  Due  to  the  perspective,  it  is  difficult  to  determine  the 

142  parameters  of  these  waves,  i.e.  their  wavelength  seems  to  decrease  when  they  are  located 

143  further  from  the  observer.  The  final  processed  images  will  permit  to  make  these 


144  measurements  using  any  well-known  eomputational  teehnique,  such  as  a  Fast  Fourier 

145  Transform. 

146 

147  3.  Processing  technique 

148  Before  processing  an  NLC  image,  it  is  necessary  to  know  the  observation  parameters 

149  of  the  imaging  system.  The  determination  of  the  horizontal  and  vertical  fields  of  view 

150  (FOVx  and  FOVy)  and  the  azimuth  acam,  elevation  Pcam  and  rotation  Rotcam  of  the  optical 

151  axis  of  the  camera  is  done  with  a  well-proven  star  calibration  method  (Garcia  et  ah,  1997). 

152  Knowing  the  exact  geographical  position  of  the  observer  and  the  exact  time  of 

153  observation,  it  is  possible  to  identify  the  stars  visible  in  a  “calibration”  image  acquired 

154  when  the  sky  is  dark  enough,  usually  at  the  beginning  or  at  the  end  of  the  NLC  season. 

155  The  imager  parameters  can  then  be  extracted  by  applying  a  least-square  algorithm  as  the 

156  exact  star  locations  are  known. 

157  Due  to  the  observation  geometry,  the  shape  of  a  projected  square  pixel  won’t  be 

158  another  square  pixel  and  its  size  will  vary,  depending  on  its  elevation  and  azimuth. 

159  Therefore,  it  is  necessary  to  process  each  image  backward,  starting  from  a  pixel  of  the 

160  final  image  and  looking  for  the  corresponding  intensity  in  the  original  data  (Garcia  et  ah, 

161  1997;  Pautet  and  Moreels,  2002). 

162  Subsequently,  we  need  to  determine  the  parameters  of  the  final  projected  image, 

163  starting  with  its  size.  It  is  given  by  the  difference  between  the  ranges  corresponding  to  the 

164  highest  and  lowest  visible  points  of  the  NLC  layer.  The  elevation  ^Max  of  the  top  of  the 

165  image,  given  by 


PMax=Pca»+FOVfy2 


16^1) 


167  Correspond  to  the  elosest  visible  point,  loeated  at  a  distanee  Dmin  from  the 

168  observation  site 

^Min  ~(^E  ^^NLc))  1691^2) 

170  Where  RMin  is  defined  as 

Rm™  =  ('R-E-tanPMax  +  (tan^PMax-CR-E  +hNLc)^  +hNLC-(2-RE  +hNLc))’^^)/(l  +  tan^PMax)^ 

171  (3) 

172  With  Re  the  average  Earth  radius,  taken  equal  to  6371  km, 

173  and  hNLC,  the  altitude  of  the  NEC  layer,  taken  equal  to  82.5  km  in  this  analysis. 

174  The  elevation  PMin  corresponding  to  the  furthermost  point  has  to  be  ehosen  depending 

175  on  the  display  itself.  In  many  instanees,  wave  struetures  are  not  resolved  in  the  lower  part 

176  of  the  images  due  to  either  the  inereased  sky  brightness  or  the  deereased  resolution,  henee, 

177  ehoosing  the  bottom  image  pixel  as  the  furthermost  point  is  not  always  relevant. 

178  Choosing  a  higher  minimum  elevation  also  limits  the  size  of  the  projeeted  image  and 

179  improves  the  resolution  of  the  wave-like  struetures  in  the  final  proeessed  produet. 

180  The  maximum  distanee  Dmax  is  then  ealeulated  using  the  same  equation  as  for  Dmin 

Dwax  =  (Re  +  hNLc)-asin(RMax/(RE  +  ^nlc))  18^) 

1 82  With  RMax  defined  as 

Rfviax  =(-RE-tanPMin  +  (tan ^ P .(R ^  +hj^Lc)^  +hNLc-(2-RE  +  hNLc))’^^)/(l  +  tan^  P^in)^ 

183  (5) 

184  The  next  step  is  to  determine  the  spatial  resolution  of  the  final  image  in  order  to  fit 

185  the  entire  projeeted  area.  The  vertieal  size  (in  meters)  depends  on  the  distanee  between 

186  the  minimum  and  the  maximum  ranges,  DMin  and  Djviax,  ealeulated  in  (2)  and  (4).  The 


1 87  final  value  Disty  is  given  by 

188  Disty  =100.[int((D^^^  - D^^J/100)  + 1] 


(6) 
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Where  int(x)  represents  the  integer  part  of  the  number  x. 

The  resolution  (in  meters  per  pixel)  is  obtained  by  dividing  the  vertieal  size  of  the 
image  Disty  by  the  number  of  pixels  Yimg  of  the  final  image 

Seale  =  DistY/Yj^g  (7) 

The  pixel  coordinates  (Xq,  Yq)  of  the  observation  site  in  the  projected  image  are 

Xo=X,„^/2  (8) 

Yo  =((D^,,  -D^J/2-Scale.Y,,g  /2)/Scale  +  Y,„^  (9) 

Here  Ximg  and  Yimg  represent  the  horizontal  and  vertical  pixel  sizes  of  the  image,  with 
the  origin  point  (0,  0)  corresponding  to  the  top  left  pixel. 

The  distance  Disto  (in  meters)  between  this  origin  point  and  the  border  of  the  image  is 
given  by 

Dist  .=  ((D  „„+  D  Scale  -  Y„  >2  (10) 

If  we  consider  a  point  A  in  the  NLC  layer,  its  coordinates  on  the  projected  image  will  be 
noted  (xa,  yA),  in  the  pixels  scale  and  (Xa,  Ya)  in  the  metric  scale,  with  the  observation 
site  as  the  origin.  The  metric  scale  coordinates  of  this  point  are  given  by 

X,=(x,-X,„,)'2.Scale  (11) 

Y^=(Dist,+Y,„,-yJScale  (12) 

The  azimuth  angle  Ua  between  the  camera  axis  and  the  direction  of  A  (see  Fig.  2)  is 

=atan(XAAfA)  (13) 

Garcia  et  al.  (1997)  have  already  calculated  the  elevation  angle  Pa  of  A  above  the  horizon 
for  the  similar  case  of  airglow  observations.  Fig.  3  is  a  modification  of  their  Fig.  6, 
adapted  to  the  current  situation. 


211 

The  distance  r  between  the  projection  of  the  observation  site  0  and  the  point  A  at  the 

212 

NLC  altitude  is 

213 

■■=A*’+y2 

(14) 

214 

The  other  parameters  shown  in  Fig.  3  are  given  by 

215 

(p  —  r/(R  E  +  h  ^Ec  ) 

(15) 

216 

d  =  2(RE  +hNEc)sin^/2 

(16) 

217 

(j)  =  {Tl-  (p)l2 

(17) 

218 

l  =  (hNLc  +d"  -2h^Ecdcos^) 

(18) 

219 

Finally  the  elevation  Pa  is  equal  to 

220 

221 

J dsin^  ^ 

Pa  =  cos 

V  1  y 

(19) 

222 

For  the  low  elevations,  the  atmospheric  refraction  has  to  be  taken  into  account.  In  our 

223 

case,  the  Lagrange  equation  is  used  for  an  elevation  higher  than  13°.  The  refraction  angle 

224 

is  then  given  by 

R  =  «o(l-Po)tanpA  -«o 

r  a  ^ 

^  2  ; 

tan'pA 

(18) 

225 

For  normal  temperature  and  pressure  conditions 

226 

ao  =  0.00029255 

227 

Po  =  0.001254  then  a  = 

Po-ao/2  = 

0.001108 

228 

At  lower  elevations,  the  atmospheric  density  distribution  is  more  complex  and  the 

229 

refraction  angle  is  approximated  by  the  following  equation 

R  =  «o  1-^  A  -^inpA^ 
1  2  J\a 

erf(X-ln/7j 
2  7 

(19) 
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With  the  function  erf(X)  equal  to 

(20) 
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233 
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240 

241 

242 

243 

244 


X  00  'Y'2k+1 

erf(x)  =  f  e-“'  du  =  - 

6^k!(2k  + 


+  1) 


And 


X  =  cosP^  / 42a 


(21) 


Hnlc 

A=e“'’  (22) 

Hp=  10400  m 

Fig.  4  shows  the  point  A’,  projection  of  A  on  the  CCD  plane  (the  original  picture).  Its 


coordinates  (Xa’,  Ya’),  in  meters,  are  given  by 


Y  cospc,^tanp^  -sinp^am  cosa^ 
cospcamcosa^ +sinpcamtana^ 


X 


A 


ftana^cospcam 


f  cospcamtanpA  -  sin  P Cam  COSOa 
cospca^cosa^  +sinpc,„tanaA 


sinPcam 


Where  f  is  the  focal  length  of  the  camera  lens. 


(23) 

(24) 


Finally,  the  pixel  coordinates  (xa’,  yA’)  of  the  point  A’  on  the  CCD  plane  are 


XA’  =  X^'.XScale.cos(-Rotc^^)  +  YA’.YScale.sin(-Rotc,„)  +  X;^g/2  (25) 

y A  -  X A  '-Xs^^i^  .sin(-Rotc,„ )  -  Y^  '.Ys^^,,  .cos(-Rotc,„  )  +  Yi^g/2  (26) 

With  Xscaie  and  Yscaie  the  horizontal  and  vertical  scales,  in  pixels  per  meters,  of  the 
original  image 

Xs„„  =  (X,„  /  2>(f.tan(FOV,,/2))  (27) 

Ys„„  =  (y„  /  2>(f.tan(FOV,,/2))  (28) 


It  is  not  necessary  to  know  the  focal  length  of  the  lens  because  the  coordinates  (xa’,  yA’) 
will  independent  of  this  parameter  which  disappears  when  integrating  equations  (23)  to 
(28)  in  a  computer  algorithm. 


245  The  final  results  are  proeessed  as  real  numbers,  so,  for  a  more  aeeurate  measurement 

246  of  the  value  of  the  pixel  A,  the  intensity  has  to  be  interpolated  using  the  pixels 

247  surrounding  the  point  A’  in  the  original  image. 

248  The  main  souree  of  error  in  the  pixel  eoordinates  is  due  to  the  uneertainty  related  to 

249  the  altitude  of  the  NLC  layer.  The  average  value  is  equal  to  82.5  km  but  it  may  vary  from 

250  80  to  85  km  (Gadsden  and  Sehroder,  1989).  If  the  altitude  differenee  is  ±2  km,  the  error 

25 1  on  the  distanee  between  the  observation  site  and  a  point  of  the  NLC  layer  will  be  ~2%. 

252  Noetilueent  elouds  are  usually  visible  when  they  are  loeated  between  100  and  1000  km 

253  away  from  the  observer,  depending  on  the  latitude  of  the  observation  site  and  the  time  of 

254  the  year.  For  these  minimum  and  maximum  ranges,  the  error  on  the  position  of  a  point  of 

255  the  NLC  layer  will  vary  between  2.4  and  13  km  (2.4%  and  1.3%),  respeetively. 

256  As  an  example.  Fig.  5  shows  the  projeeted  image  of  the  NLC  pieture  presented  in  Fig. 

257  1,  assuming  a  layer  altitude  of  82.5  km.  The  horizontal  and  vertieal  fields  of  view  are  53° 

258  and  41°,  respeetively.  The  NLC  layer  is  presented  as  seen  from  the  top,  with  respeet  to 

259  the  geographieal  orientation,  and  its  seale  is  visible  in  the  bottom  right  eomer.  The  short- 

260  period  gravity  waves  now  appear  with  their  real  horizontal  wavelengths  and  direetions  of 

261  propagation.  Their  horizontal  observed  phase  speed  may  be  measured  when  looking 

262  through  a  sequenee  of  images  and  the  intrinsie  horizontal  phase  speed  dedueed  knowing 

263  the  veloeity  and  direetion  of  the  baekground  wind.  Well-developed  tools  employing 

264  Fourier  transform  teehniques  were  used  to  measure  these  parameters  (Gareia  et  ah,  1997; 

265  Coble  et  ah,  1998).  The  results  for  the  larger-seale  wave  (denoted  by  A)  moving 

266  northeastward,  as  indieated  by  the  white  arrow,  in  the  right  part  of  the  image  are; 

267  horizontal  wavelength  kh  =  29.3±3  km,  observed  horizontal  phase  speed  Uh  =  18.3±5  m/s, 

268  observed  period  x  =  26.7±2  min,  and  direetion  of  propagation  due  to  North  0  =  31±3°.  It 


269 


is  important  to  notice  that  due  to  the  size  of  the  field  of  view  (53°x41°),  the  maximum 

270  observable  wavelength  is  ~  100km,  even  though  longer  waves  have  been  reported  at  the 

271  NLC  altitude  (Carbary  et  ah,  2001;  Chandran  et  ah,  2010). 

272  It  is  also  possible  to  measure  the  eharaeteristies  of  the  ripple  event  (noted  B),  visible 

273  on  the  eenter  right  part  of  the  projeetion:  horizontal  wavelength  kh  =  5.2±3  km,  observed 

274  horizontal  phase  speed  Uh  =  68.6±5  m/s,  observed  period  x  =  1.3±1  min,  and  direetion  of 

275  propagation  due  to  North  0  =  72±3“.  These  values  are  in  the  range  of  previous  studies  of 

276  mesospherie  ripple-like  waves  (e.g.,  Taylor  et  ah,  1997;  Heeht,  2004). 

277 

278  4,  Short-period  gravity  waves  analysis 

279  During  the  past  deeades,  atmospherie  gravity  waves  have  been  intensely  studied  using 

280  optieal  instrumentation  deteeting  their  signature  in  the  upper  atmosphere  nighttime 

281  airglow  emissions  (altitudes  '-'80-100  km).  Most  of  the  observations  were  eondueted 

282  from  sites  loeated  at  low  (e.g.,  Taylor  et  ah,  1997;  Medeiros  et  ah,  2003;  Nakamura  et  ah, 

283  2003;  Pautet  et  ah,  2005;  Wrasse  et  ah,  2006)  or  mid-latitudes  (e.g.,  Wu  and  Killeen, 

284  1996;  Heeht  et  ah,  2001;  Ejiri  et  ah,  2003).  The  literature  about  ground-based  short- 

285  period  gravity  wave  optieal  studies  at  high-latitudes  is  not  very  extensive  (Clairemidi  et 

286  ah,  1985;  Taylor  and  Henriksen,  1989;  Dalin  et  ah,  2004;  Nielsen  et  ah,  2006;  Nielsen  et 

287  ah,  2009;  Bageston  et  ah,  2009)  due  to  the  diffieult  observation  eonditions  like  the 

288  presenee  of  twilight  during  the  summer  months  and/or  the  possible  auroral  eontamination 

289  of  the  faint  optieal  airglow  emissions.  For  these  reasons,  the  data  deseribed  herein 

290  provide  a  broad  and  unique  dataset  to  investigate  the  small-seale  gravity  waves 

291  propagating  through  the  MLT  region  at  a  high-latitude  site  (59.4°N)  during  the  summer 


292 


months. 


293  Between  2004  and  2008,  this  imager  has  run  automatieally  almost  every  night  from 

294  May  to  September,  eovering  the  whole  loeal  NLC  season.  If  we  eonsider  the  period 

295  extending  from  the  day  of  the  earliest  sighting  to  the  day  of  the  latest  sighting  during 

296  these  5  years  (May  31-June  1  to  August  14-15),  the  eamera  has  operated  during  322 

297  nights.  247  of  these  nights  (77%)  were  not  totally  eloudy,  allowing  seeing  possible  NLC 

298  displays,  and  during  107  of  them  (43%  of  the  observable  nights)  noetilueent  elouds  were 

299  aetually  visible.  From  these  107  nights,  32  have  been  analyzed  to  determine  the 

300  eharaeteristies  of  the  gravity  wave  events  appearing  in  the  images.  During  the  other 

301  nights,  the  eloud  eoverage  was  too  extensive  to  deteet  any  NLCs  or  the  NLCs  were  too 

302  weak  or  too  far  away  from  the  observation  site  to  resolve  any  wave  strueture. 

303  After  proeessing  the  images  as  deseribed  in  the  first  part  of  this  paper,  the  wave-like 

304  patterns  were  isolated  and  their  parameters  (direetion  of  propagation,  horizontal 

305  wavelength,  observed  horizontal  phase  speed  and  observed  period)  proeessed  using  a 

306  well-proven  Fast  Fourier  Transform  teehnique  (Gareia  et  ah,  1997;  Coble  et  ah,  1998). 

307  These  wave  events  were  divided  in  two  groups,  depending  on  their  aspeet  and 

308  eharaeteristies:  the  shortest-wavelength,  loealized  and  short-lived  waves  (37  events), 

309  often  referred  as  “ripples”  or  “billows”,  and  the  longer- lived,  more  spread  waves 

310  exhibiting  longer  duration  and  wavelengths  (>10  km)  referred  as  “bands”  (30  events). 

311  The  ripple  waves  are  not  deseribed  any  further  in  this  study  beeause  they  are  mostly  due 

312  to  loealized  instabilities  (e.g.,  Taylor  and  Hapgood,  1990;  Fritts  et  ah,  1993,  1994;  Heeht 

313  et  ah,  1997)  and  are  not  thought  to  provide  direet  information  on  freely  propagating 

314  gravity  waves  (Heeht,  2004).  The  seeond  group  of  events  eorresponds  to  the  short-period 

315  (<  Ih)  gravity  waves  propagating  from  the  troposphere  through  the  MLT  region  and  that 


316  are  believed  to  be  responsible  for  the  momentum  deposition  and  part  of  the  global 

317  eireulation  in  the  upper  atmosphere  (e.g.,  Fritts  and  Alexander,  2003). 

318  The  horizontal  wavelength  (Fig.  6a),  the  observed  horizontal  phase  speed  Ux 

319  (Fig.  6b)  and  the  resulting  observed  period  x  (Fig.  6o)  are  showed  in  histogram  plots  to 

320  summarize  the  charaeteristies  of  the  30  short-period  gravity  wave  events  observed  in  the 

321  NLC  images  between  2004  and  2008.  The  mean  values  are  respeetively  25.1±1  km,  with 

322  60%  of  the  wavelengths  between  20  and  30  km,  27.2±5  m/s  with  83%  of  the  phase 

323  speeds  between  10  and  40  m/s,  and  21.9±2  min  with  60%  of  the  periods  between  10  and 

324  30  min.  The  eomparison  with  previous  long-term  studies  shows  that  the  average 

325  horizontal  wavelength  is  very  similar  for  gravity  waves  observed  in  the  airglow  emissions, 

326  at  slightly  higher  altitudes,  than  in  the  NC  layer.  For  example,  Taylor  et  al.  (1997)  (2°S), 

327  Nakamura  et  al.  (2003)  (7°S)  and  Medeiros  et  al.  (2003)  (23°S),  measured  values  of  24 

328  km,  29  km  and  23  km,  respeetively,  for  sites  loeated  at  low  latitudes.  At  mid-latitude 

329  loeations,  the  results  are  similar  with  horizontal  wavelengths  reported  by  Heeht  et  al. 

330  (2001)  (40°N)  and  Ejiri  et  al.  (2003)  (44°N)  of  25  and  24  km,  respeetively.  In  a  long-term 

331  study  at  a  high-latitude  site,  condueted  during  the  2000  and  2001  austral  winters  from 

332  Halley,  Antaretica  (76°S),  Nielsen  et  al.  (2009)  obtained  an  average  horizontal 

333  wavelength  of  26±lkm.  During  the  2007  austral  winter  season  (April  to  October  2007), 

334  Bageston  et  al.  (2009)  observed  an  horizontal  wavelength  varying  from  10  to  60  km,  with 

335  a  maximum  occurrence  between  15  and  35km,  at  Comandante  Ferraz,  Antarctica  (62°S). 

336  All  those  results  are  in  very  good  agreement  and  confirm  the  systematic  global  presence 

337  of  short-period  gravity  waves  at  mesospheric  altitude.  In  contrary,  the  observed 

338  horizontal  phase  speeds  obtained  from  Sweden  are  different  from  the  previous  studies 

339  with  much  smaller  values  (average  27.2±5  m/s).  In  the  past  measurements,  the  average 


340  phase  speeds  vary  signifieantly  from  one  site  to  another,  without  any  apparent  relation 

341  with  the  latitude.  At  low-latitudes,  Taylor  et  al.  (1997)  (2°S)  found  a  mean  value  of  48 

342  m/s  while  Nakamura  et  al.  (2003)  (7°S)  obtained  56  m/s  and  Medeiros  et  al.  (2003) 

343  (23°S)  only  26  m/s.  The  numbers  can  also  diverge  greatly  at  mid-latitudes;  for  example 

344  Hecht  et  al.  (2001)  (40°N)  measured  an  average  value  of  50  m/s  and  Ejiri  et  al.  (2003) 

345  (44°N)  only  35  m/s.  The  comparison  between  the  Swedish  dataset  and  the  measurements 

346  made  at  Halley,  Antarctica  (76°S),  is  even  more  striking  with  an  average  value  of  48  m/s 

347  for  the  Southern  location,  almost  twice  the  value  obtained  in  the  current  study.  However, 

348  Clairemidi  et  al.  (1985)  who  operated  a  near-IR  photographic  camera  at  Sodankyla, 

349  Finland  (67.3°N),  reported  observed  phase  speeds  of  15.6  m/s  varying  from  7  to  26  m/s, 

350  during  a  limited  10-day  study  in  Dec  1981 -Jan  1982.  These  results  have  to  be  taken 

351  carefully  because  comparison  between  the  datasets  would  only  be  possible  using  the 

352  intrinsic  horizontal  phase  speeds,  but  in  most  cases,  this  parameter  has  not  been 

353  processed  due  to  a  lack  of  knowledge  on  the  wind  velocity  at  the  observation  site. 

354  Nevertheless,  this  comparison  gives  a  good  estimate  of  the  differences  between  the 

355  studies,  especially  for  similar  latitudes  and  period  of  the  year.  The  observed  periods, 

356  which  are  directly  related  to  the  horizontal  phase  speeds,  are  also  spread  over  a  broad 

357  range  of  values  (from  8  to  15  min)  and  are  smaller  in  the  previous  studies  than  for  the 

358  current  measurements  (21 .9±2  min). 

359  Fig.  7  shows  the  azimuthal  distribution  of  the  observed  waves.  There  exists  a 

360  strong  anisotropy  with  a  majority  of  the  waves  propagating  towards  the  north  or  the 

361  north-east  (60%  of  the  headings  between  330°  and  60°).  Such  asymmetry  has  been 

362  observed  before,  from  both  hemispheres,  with  a  preferential  poleward  motion  during  the 

363  summer  months,  usually  slightly  tilted  towards  the  east  (e.g.,  Taylor  et  ah,  1997;  Hecht  et 


364  al,  2001;  Ejiri  et  al,  2003;  Medeiros  et  al,  2003).  Two  main  reasons  have  been  proposed 

365  to  explain  this  anisotropy:  Taylor  et  al.  (1993)  demonstrated  that  the  baekground  wind 

366  can  alter  the  upward  flux  of  waves  depending  upon  their  phase  speeds  and  headings.  If 

367  the  wind  speed  component  parallel  to  the  direction  of  propagation  of  a  wave  becomes 

368  equal  to  the  observed  phase  speed  of  this  wave,  the  wave  can  be  absorbed  into  the 

369  background  medium  and  consequently  prohibited  in  reaching  higher  altitudes.  This  so- 

370  called  critical  layer  filtering  mechanism  may  explain  the  strong  anisotropy  in  the  wave 

371  propagation  seen  from  Stockholm.  Localized  wind  measurements  throughout  the  lower 

372  and  middle  atmosphere  are  rare  and  do  not  exist  for  the  time  and  location  considered  here. 

373  To  address  the  potential  impact  of  the  intervening  wind  field,  we  utilize  derived  winds 

374  from  an  advanced  numerical  weather  prediction  system,  the  Navy  Operational  Global 

375  Atmospheric  Prediction  System  -  Advanced  Level  Physics  High-Altitude  (NOGAPS- 

376  ALPHA)  forecast  model  and  data  assimilation  system,  which  has  recently  been  extended 

377  into  the  lower  MLT  region  (Hoppel  et  ah,  2009;  Lckermann  et  ah,  2009).  Localized 

378  NOGAPS-ALPHA  winds  from  69°N  has  been  validated  against  meteor  radar  winds 

379  (Stevens  et  ah,  2009).  In  this  paper  we  use  winds  from  63°N  and  18°E,  which 

380  corresponds  to  the  region  north  of  the  observation  where  majority  of  the  NLCs  were 

381  observed. 

382  Lig.  8  shows  the  meridional  and  zonal  background  wind  components  for  July 

383  2007  at  0:00  UT.  The  strong  westward  zonal  wind  just  below  the  NLC  layer  (60  m/s) 

384  may  justify  the  small  eastward  component  in  the  gravity  waves  propagation  but  critical 

385  layer  filtering  due  to  the  limited  meridional  wind  (<10  ms  below  NLC  altitude)  cannot 

386  account  for  the  strong  northward  heading  although  it  may  explain  the  fact  that  so  many 

387  small  phase  speed  events  are  not  filtered  out. 


388  Another  possible  reason,  and  more  likely  reason,  for  the  headings  anisotropy  is 

389  the  non  uniformity  in  the  distribution  of  the  souree  regions.  Atmospheric  gravity  waves 

390  have  different  origins:  orographic  forcing  by  mountains,  frontal  zones,  or  tropospheric 

391  convection.  During  the  summer  months,  localized  convection  develops  south  of  Sweden, 

392  mainly  over  central  and  western  European  countries;  in  addition,  the  regions  further  north 

393  exhibit  weak  or  no  convection.  Short-period  gravity  waves  originating  from  this  type  of 

394  sources,  which  have  been  previously  observed  in  the  mesospheric  airglow  over  Europe 

395  (Taylor  and  Hapgood,  1988),  can  be  accounted  for  the  anisotropy  in  the  propagation 

396  direction.  Moreover,  Holton  et  al.  (2002)  showed  that  the  dominant  vertical  wavelength 

397  of  the  short-period  gravity  waves  generated  by  convection  depends  on  the  characteristics 

398  of  the  source  region.  The  observed  wave  phase  speeds  are  also  expected  to  be 

399  proportional  to  the  vertical  wavelength  (through  the  gravity  wave  dispersion  relationship) 

400  and  hence  to  the  heating  depth.  Thus,  the  limited  heating  at  these  latitudes  associated  with 

401  less  intense  convection  compare  to  observation  sites  located  at  lower  latitudes,  could 

402  have  excited  waves  with  smaller  horizontal  phase  speeds.  In  order  to  make  sure  that  the 

403  gravity  waves  observed  from  Sweden  are  generated  by  close  tropospheric  convection,  it 

404  would  be  necessary  to  determine  the  nature  of  the  wave  propagation  (freely  propagating 

405  or  ducted  waves)  and  associate  each  wave  with  a  potential  source  region  through  a  ray- 

406  tracing  technique. 

407 

408  5,  Summary 

409  Noctilucent  clouds  have  been  observed  and  studied  for  more  than  100  years.  Their 

410  exact  composition  and  the  conditions  of  their  formation  are  still  a  subject  of  ongoing 

411  research;  nevertheless,  they  are  an  interesting  tool  in  order  to  investigate  the  atmospheric 


412  dynamics,  especially  the  short-period  gravity  waves  propagating  through  the  MLT  layer. 

413  Standard  nighttime  airglow  measurements  are  not  possible  at  high-latitude  during  the 

414  summer  months  beeause  of  the  persistent  twilight  and  diffieult  during  the  winter  months 

415  due  to  the  auroral  eontamination.  The  observation  of  gravity  waves  using  noetilueent 

416  eloud  images  can  be  seen  as  a  eomplementary  method  for  these  latitudes. 

417  Sinee  2004,  a  eamera  running  automatieally  every  summer  night  on  the  Stoekholm 

418  University  eampus  has  obtained  a  large  dataset  of  NLC  pietures,  some  of  them  exhibiting 

419  gravity  wave  signatures.  The  teehnique  deseribed  in  this  paper  was  developed  to  proeess 

420  these  NLC  images  in  order  to  easily  measure  the  parameters  of  the  wave  events  they 

421  display.  This  dataset  eonfirm  the  presenee  of  short-period  gravity  waves  at  high-latitude 

422  sites  during  the  summer  months  and  the  possibility  to  suoeessfully  study  their 

423  eharaeteristies,  under  otherwise  diffieult  eonditions.  The  parameters  of  30  wave  events 

424  have  been  measured  and  eompared  with  preeeding  results  obtained  using  airglow 

425  observations,  mainly  from  low  or  mid-latitude  sites.  The  distribution  of  the  horizontal 

426  wavelengths  is  similar  to  previous  studies  but  the  observed  horizontal  phase  speeds  are 

427  slower  than  for  the  other  measurements.  This  faet  eannot  be  aeeounted  for  through 

428  eritieal  layer  filtering  but  may  instead  be  a  result  of  weaker  souree  foreing  as  eompared  to 

429  low-  and  mid-latitudes.  Furthermore,  the  direction  of  propagation  of  the  wave  events 

430  exhibits  a  strong  anisotropy  with  a  majority  of  the  waves  heading  towards  the  NNE.  This 

431  poleward  tendeney  during  the  summer  months  has  been  previously  notieed  from 

432  numerous  sites.  In  the  ease  of  the  gravity  waves  observed  above  Sweden,  the  distribution 

433  of  the  souree  regions  might  be  a  more  preponderant  factor  than  the  filtering  due  to  the 

434  baekground  wind  field.  The  limited  period  of  observation  eaeh  year  (only  the  NLC 

435  season  from  June  to  August)  does  not  allow  following  the  progression  of  the  wave 


436  headings  or  the  observed  horizontal  phase  speeds  as  a  funetion  of  the  variations  in  the 

437  baekground  wind  or  the  evolution  of  the  gravity  wave  sourees.  The  souree  region  of  the 

438  waves  is  eurrently  under  investigation  using  ray-traeing  teehnique.  The  results  should 

439  permit  to  determine  the  prevailing  nature  of  the  gravity  wave  motions.  Further 

440  eomparison  between  eurrent  and  additional  wave  measurements  and  wind  data  are  also 

441  needed  to  eonfirm  the  eurrent  study  and  to  elarify  the  role  of  the  gravity  waves  in  the 

442  dynamies,  the  momentum  deposition  and  the  global  eirculation  in  the  upper  atmosphere 

443  at  high-latitudes,  and  their  effeets  on  noetilueent  elouds  generation  and  fading.  Finally,  it 

444  would  be  of  great  interest  to  eompare  the  gravity  waves  observed  in  NLC  images  from 

445  the  ground  with  the  data  obtained  by  satellite-based  instruments  like  CIPS  onboard  the 

446  AIM  satellite  (Chandran  et  ah,  2009)  and  with  other  high-latitude  datasets. 
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599  Figure  captions 

600  Fig.  1.  Noctilucent  cloud  image  taken  from  Frescati,  Sweden  (59.4°N),  on  July  16^'^,  2005, 

601  at  22:30  UT  (0:30  LT). 

602  Fig.  2.  Determination  of  the  azimuth  angle  Ua  of  a  point  A  of  the  NLC. 

603  Fig.  3.  Plot  showing  the  NLC  layer  geometry. 

604  Fig.  4.  Geometry  of  the  projection  of  a  point  of  the  NLC  layer  on  the  CCD  plane. 

605  Fig.  5.  Projection  of  the  NLC  image  presented  in  Fig.  1.  Note  the  short-period  gravity 

606  wave  A  on  the  right  part  of  the  image  and  the  ripples  B  close  to  the  right  edge. 

607  Fig.  6:  Histogram  plots  showing  the  distribution  of  horizontal  wavelengths  (a),  observed 

608  horizontal  phase  speeds  (b)  and  observed  periods  (c)  of  the  short-period  gravity  wave 

609  events  observed  between  2004  and  2008  from  Stockholm,  Sweden 

610  Fig.  7:  Azimuthal  plot  showing  the  propagation  direction  of  the  waves  for  the  summers 

611  2004  to  2008.  Note  the  anisotropy  with  the  preferential  northward  heading. 

612  Fig.  8:  Plot  of  the  zonal  and  meridional  wind  components  derived  from  the  NOGAPS- 

613  ALPHA  assimilation  model  for  July  2007  at  0:00UT  and  63°N. 
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Fig.  1.  Noctilucent  cloud  image  taken  from  Freseati,  Sweden  (59.4°N), 
on  July  16^  2005,  at  22:30  UT  (0:30  LT). 
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639  Fig.  2.  Determination  of  the  azimuth  angle  Ua  of  a  point  A  of  the  NLC. 
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658  Fig.  4.  Geometry  of  the  projeetion  of  a  point  of  the  NLC  layer  on  the  CCD  plane. 
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Fig.  5.  Projection  of  the  NLC  image  presented  in  Fig.  1.  Note  the  short-period  gravity 
wave  A  on  the  right  part  of  the  image  and  the  ripples  B  close  to  the  right  edge. 
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Fig.  6;  Histogram  plots  showing  the  distribution  of  horizontal  wavelengths  (a),  observed 


horizontal  phase  speeds  (b)  and  observed  periods  (c)  of  the  short-period  gravity  wave 


events  observed  between  2004  and  2008  from  Stockholm,  Sweden 
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Fig.  7;  Azimuthal  plot  showing  the  propagation  direction  of  the  waves  for  the  summers 
2004  to  2008.  Note  the  anisotropy  with  the  preferential  northward  heading. 
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717  Fig.  8:  Plot  of  the  zonal  and  meridional  wind  components  derived  from  the  NOGAPS- 

718  ALPHA  assimilation  model  for  July  2007  at  0;00UT  and  63°N. 
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